Diffusion of single ellipsoids under quasi-2D confinements 
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We report video-microscopy measurements of the translational and rotational Brownian motions 
of isolated ellipsoidal particles in quasi-two-dimensional sample cells of increasing thickness. The 
long-time diffusion coefficients were measured along the long (Da) and short (Db) ellipsoid axes, 
respectively, and the ratio, Da/Dt,, was determined as a function of wall confinement and particle 
aspect ratio. In three-dimensions this ratio (Da/Dt) cannot be larger than two, but wall confinement 
was found to substantially alter diffusion anisotropy and substantially slow particle diffusion along 
the short axis. 
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INTRODUCTION 

In many biological and industrial processes, diffusing 
particles are non-spherical and move in confined geome- 
tries. Examples of particles in this scenario include pro- 
teins diffusing in membranes [Ij and very fine grains mi- 
grating through pores in porous media. To date, quan- 
titative measurements of anisotropic particle diffusion in 
confined geometries have been limited. However, new 
particle fabrication and imaging technologies combined 
with new image analysis tools now make the direct mea- 
surement of the diffusion of anisotropic particles readily 
possible. Thus, in this contribution we investigate the 
anisotropic diffusion of isolated ellipsoidal particles con- 
fined between two parallel plates. 

The Brownian diffusion coefficient D of an isolated 
spherical particle is well understood. It is inversely pro- 
portional to the drag (or friction) coefficient 7 via the 
Einstein relation. 



D = keTh 



(1) 



where fc^ is the Boltzmann constant and T is the temper- 
ature. For a prolate spheroid with long axis of length 2a 
and two short axes of length 26, translational diffusion 
is anisotropic and is described by diffusion coefficients 
Da = ksT/^a along the long axis, and Db — ksTl^b 
along the short axes. The rotational diffusion coeffi- 
cient of the prolate spheroid about its short axes is 
De = kBT/^g. Generally, the drag coefficients 7a, 7fc 
and 70 depend on the shape and size of the ellipsoid. 
Brownian motion of anisotropic particles was first seri- 
ously considered by F. Perrin 0, y] who computed these 
drag coefficients analytically for a spheroid diffusing in 
three dimensions (3D). Interestingly, the ratio Da/Db 
varies from one to two in 3D, as the spheroid aspect ratio 
4> = a/h varies from one to infinity. 

The problem of diffusion in confined geometries, such 
as quasi-2D media, is different from the 3D case as a re- 
sult of a complex interplay between hydrodynamic drag. 



the boundaries of the medium, and the particle geom- 
etry. Surfaces near a moving particle modify fluid flow 
fields, often increasing particle hydrodynamic drag. A 
full theoretical formulation of wall hydrodynamic effects 
has been developed for one sphere (or ellipsoid) coupled 
to one wall [4]. However, for more complicated situa- 
tions, such as a sphere or an ellipsoid confined by two 
parallel walls, the only available analytical solutions are 
for weak confinement in a few special symmetric config- 
urations Recent numerical calculations Q, on the 
other hand, have been developed to derive the hydrody- 
namic drag of a single sphere and a linear chain of spheres 
confined more strongly in quasi-2D. 

On the experimental side, the hydrodynamic drag of 
single spheres in weak confinement have been measured 
@, 13 J and video microscopy has been applied recently 
to measure anisotropic particle diffusion, including ellip- 
soids in quasi-2D [1] and 3D 0, colloidal clusters near 
one wall [lO| , and carbon nanotubes in weak conflnement 
111 . In the present contribution we report measurements 
of hydrodynamic drag on ellipsoids in quasi-2D, conflned 
between two parallel walls. We explore the strong con- 
finement regime where drag coefficients are not readily 
available from theory and simulation, and we report on a 
light interference method to accurately measure the con- 
finement. We find that the diffusion anisotropy is made 
stronger and the diffusion along ellipsoid short axes is 
dramatically slowed due to wall confinement. The exper- 
iment and analyses are similar to a previous paper 
However the scope of the present work is different, focus- 
ing instead on how confinement affects diffusion, rather 
than on the detailed time-dependent Brownian dynamics 
of a single ellipsoid with the greatest diffusion anisotropy. 



THEORY BACKGROUND 

When a spheroid with semi-axes (a, h, b) moves along 
one of its principle axes with velocity v, through an un- 
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FIG. 1: The geometric factors G in Eq. (|3]) as a function of 
aspect ratio (j)- 



bounded quiescent fluid with viscosity ry at low Reynolds 
number, then the translational and rotational (about 
short axis) drag coefhcients affecting the spheroid are 



7 — 6Trr]bG, 
76) = driVGe- 



(2a) 
(2b) 



V is the volume of the spheroid and G is the geomet- 
ric factor that renders the ellipsoid different relative to 
the case of a sphere. The geometric factors for prolate 
spheroids diffusing in 3D are analytically given by Per- 
rin's equations 
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Here = a/6 is the aspect ratio. When 0=1, then 
G ~ Gg = 1, and Eq. ([2]) reduces to the translational 
and rotational Stokes laws for a sphere. Note also that 
Eqs. ^ and ^ are obtained using stick boundary con- 
ditions, valid when the particle is much larger than fluid 
molecules [11, El. In Fig. [2 Eq. © is plotted out as a 
function of 4> for (j> less than 10. When the aspect ratio 
</)>!, Eqs. HI), jg) and ^ yield 

fcsTln^ ^ kBTlncj) 



Da = 



2Trr]a 



(4) 



The ratio between these diffusion coefficients along long 
and short axes, i.e. Da/ Db = Gb/Ga, increases monoton- 
ically from one to two as </> increases from one to infinity 
(in 3D). In quasi-2D, however, Da/Db can be larger than 
two. 



EXPERIMENT 

The diffusion of micrometer size PMMA (polymethyl 
methacrylate) and PS (polystyrene) ellipsoids was mea- 
sured in water confined between two glass walls. Both 
PS and PMMA ellipsoids are synthesized by the method 
described in Ref. Briefly, we placed 0.5% (by 

weight) PS spheres into a 12% (by weight) aqueous PVA 
(polyvinyl alcohol) solution residing in a Petri dish. Af- 
ter water evaporation, the PVA film was stretched at 
130°C. The PS (or PMMA) spheres embedded in the 
film are readily stretched because their glass transition 
temperatures are below 130°C. After cooling to room 
temperature, the PVA was dissolved and ellipsoids ob- 
tained. Note, the initial PMMA or PS spheres must not 
be cross-linked, otherwise they cannot be stretched. We 
measured the size of ellipsoids by SEM and by optical 
microscopy. 

The ellipsoid solutions were cleaned and stabilized with 
7 mM SDS (sodium dodecyl sulfate). The ellipsoids were 
not expected to have strong interactions with the glass 
surfaces, because the solution ionic strength was more 
than 0.1 mM and the Debye screening length for the par- 
ticles was correspondingly less than 30 nm. However, it 
is difficult to estimate the ionic strength accurately in a 
thin cell because the glass surfaces can release Na+ ions 
[l^. Nevertheless, we found that the addition of 2 mM 
salt to the solution did not induce a detectable change in 
particle diffusion coefficients. This observation suggests 
that the double layers are not significantly affecting par- 
ticle diffusion. 

Glass surfaces of the sample cell were rigorously 
cleaned in a 1:4 mixture of hydrogen peroxide and sul- 
furic acid by sonication. Then the glass was thoroughly 
rinsed in deionized water and quickly dried with an air 
blow gun. Typically 0.3 fiL solution spread over the en- 
tire 1.8 X 1.8 mm^ coverslip area, and ellipsoids did not 
stick to the surfaces. Because the gravitational height, 
ksT /mg, is much larger than the cell thickness, H , the 
ellipsoids were readily suspended around mid-plane be- 
tween the two walls. Finally, the cell was sealed with UV 
cured adhesive (Norland 63). 

We measure the wall separation by light interference. 
When the cell thickness is below a few micrometers, then 
the interference colors produced by reflections from the 
two inner surfaces of the sample walls in white light illu- 
mination can be observed by eye or in the reflection mode 
of microscope, see Fig. [2]4_,B. When the wall separation 
H = 0, the effective light path difference is Al = A/2 
due to the tt phase shift of reflection at the lower sur- 
face. Thus all wavelength components of the white light 
yield a dark black color in interference at 7J = 0. When 
H > 0, the reflection light in the normal direction is a 
mixture of light with various wavelengths, and different 
wavelengths contribute with different weights to the ob- 



3 



served color. White light interference from a wedge, for 
example, will be bands of colors as in the Michel-Levy 



Chart [17[. By comparing the observed color with the 



Michel-Levy Chart, we can effectively read out the cor- 
responding Al and obtain H = Al/(2nw), where n„ is 
the refractive index of water. In the Michel-Levy Chart, 
the color starts from black at = and changes from 
red to blue periodically with period Al = 625 nm. To 
avoid misreading the color by one or more periods, we 
either made a reference wedge or we put dilute spacer 
spheres with known diameter between the glass slides 
to establish a reference thickness. Also, color bands may 
shift slightly because the illumination light is not an ideal 
white light source. This error however, should be less 
than 625/4 nm, so that the error of H is less than SH (625 
nm/4)/(2nm) = 60 nm. Although the absolute value of 
H may be subject to ~ 60 nm uncertainty as described 
above, the relative values of different H in one cell should 
be more accurate (~ 30 nm) because we can easily distin- 
guish more than 8 different colors in one band including 
deep red, light red, orange, light orange, yellow etc. 

Usually our sample thickness had less than 20 nm vari- 
ation in the central 1 mm^ area and had 1-2 variation 
over the whole 18 x 18 mm^ area. Thus, we can study the 
diffusion of ellipsoids at different H in one cell. The inter- 
ference between reflected light from top inner surface of 
the wall and the ellipsoid's top surface give rise to differ- 
ent colors (see Fig.[2|). As is the case with Newton's rings, 
the interference colors due to the two ellipsoid tips and 
the center of the ellipsoid were different. We found that 
the color only fluctuated near the two ellipsoid tips; the 
color was quite constant near the ellipsoid center. Thus 
the height fluctuation of the ellipsoids was very small and 
tumbling motions in the vertical plane were not strong. 
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FIG. 2: (color online) (A) Schematic of sample dimensions 
and the interference observation mode. (B) True interference 
color image of ellipsoid in the reflection mode of the micro- 
scope. (C) Bright field ellipsoid image in the transmission 
mode. 

Particle motions observed by microscopy were recorded 
by a CCD camera to videotape at 30 frames/sec. In the 
dilute suspension, only one ellipsoid was visible in the 
640 X 480 pixel^ = 38.4 x 51.2 ^J.m^ field of view un- 
der lOOx objective during a half-hour experiment. We 
defocus slightly so that the ellipsoid can be more accu- 
rately located along its long axis. The built-in 2D Gaus- 
sian fit function in IDL (Interactive Data Language) was 
used to locate the center and orientation of the ellipse 



in each video frame. In practice, a small percent (~3%) 
of the frames failed to be correctly tracked. Without 
these frames, the trajectory breaks into short pieces and 
very long-time behavior becomes difficult to measure. To 
capture these frames, we very slightly adjusted tracking 
parameters or image contrast and re-analyzed the im- 
ages; after these corrections roughly 3% x 3% — 0.09% 
of the frames remain incorrectly tracked. We then re- 
peated this procedure iteratively until all ~50000 frames 
in one dataset were correctly tracked. The mean square 
displacements (MSD) at time lag t = has small non- 
zero intercept due to the tracking errors. Thus we can 
estimate the spatial and angular resolution from inter- 
cepts of their corresponding MSDs p^] . The orientation 
resolution is 1°, and spatial resolutions are 0.5 pixel — 
40 nm along the particle's short axis and 0.8 pixel = 64 
nm along its long axis because of the superimposed small 
tumbling motion. 

From the image analysis, we obtained the trajec- 
tory of a particle's center-of-mass positions x(t„) = 
{x{tn),y{tn)) in the lab frame and its orientation angle 
6{tn) relative to the x-axis at times tn = n x (1/30) 
sec, see Fig. |3l We define each 1/30-sec time interval 
as a step. During the n^^ step, the particle's position 
changes by Sx{tn) = x(t„) — x(t„_i) and its angle by 
59{tn) = 6{tn) — d{tn-i)- To obtain the drag coefficients 
along long and short axes, we need to covert the mea- 
sured displacements from the fixed lab frame to the local 
body frame. Step displacements (5x„ relative to the local 
body-frame and step displacements Sxn relative to the 
fixed lab frame are related via 



COS 6n sin 9n 
— sin 9r, COS 9ri 



5x„ 
Syn 



(5) 



where 9n = (^(^n-i) +^(^n))/2, see Fig. [S] In practice, 
choosing 9n — 9{tn-i) or 9n = 9{tn) has little effect on 
our results because 9 barely changes during 1/30 s. 




FIG. 3: An ellipsoid in the x-y lab frame and the x-y body 
frame. The angle between two frames is 9{t). The displace- 
ment 5x can be decomposed as {Sx, Sy) or (Sx, 5y). 

Figure [3] shows mean-square-displacements (MSDs) of 
a 2.4 X 0.3 X 0.3 |J,m^^ ellipsoid conflned in an 846 nm 
thick cell. In both the lab and the body frame, MSDs 
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FIG. 4: (color online) Mean square displacements (MSDs) of 
a 2.4 X 0.3 X 0.3 /im^ ellipsoid confined in an 846 nm thick cell. 
Top panel: log-log plot. Bottom panel: linear plot. The four 
lines arranged from top to bottom, respectively, are MSDs 
along a, X, y and b axes. Symbols represent experimental 
data, and lines represent linear fits to these data. Insets: 
Angular MSDs. All curves exhibit diffusive behavior, and the 
diffusion coefficients, D — MSD/(2t), shown in the figure are 
derived from the best fit lines. 



are diffusive with ([Ai(t)]2) = 2Dat, {[Ay{t)]'^) =_2Dbt, 
([Ax(t)]2) = ([Ay(t)]2) = {Da + Dh)t = 2Dt and 
{{A9{t))') = 2Det. 



RESULTS AND DISCUSSION 

We repeated the experiments described above for dif- 
ferent ellipsoids under different confinement conditions. 
From the slopes of their MSDs, wc obtain Da, Di, and Dg 
of different particles as a function of confinement condi- 
tion as shown in Figs. [3 [6] and [71 respectively. Specif- 
ically, the normalized quantities, D^^ /Di = 7i/7f^, for 
i = a,b,9, are plotted as a function of increasing confine- 
ment, 2b/H. Here the 3D normalization constants D^^ 
(alternatively, 7^^), are calculated from Eqs. [2]and[3l 

Notice that 2b /H ~ corresponds to the 3D limit 
wherein D^^ /D = 1. As expected, hydrodynamic drag 
increased and the diffusion coefficients correspondingly 




2b /H 



FIG. 5: (color online) Ratio of theoretical 3D diffusion coef- 
ficient Q along the ellipsoid long axis, D^^, to the measured 
diffusion coefficient, Da, for ellipsoids confined at 2b/H. Di- 
amonds: 2.4 X 0.3 X 0.3 /xm^ (</> = 8) ellipsoids; Circles: 
3.3 X 0.635 X 0.635 /xm^ (0 = 5.2) ellipsoids; Solid circles: 
From samples with no added salt; Open circles: From samples 
with 2 mM added salt; Green solid circle: From sample with 
BSA (bovine serum albumin) covered glass surfaces; Squares: 
All other samples - lower aspect ratio spheroids with particle 
aspect ratios labeled below each data point. The accuracy 
of 2b /H for these measurements is similar to other samples. 
Dashed curves: Guides for the eye. Solid curve {(j> = 1): Re- 
plot of the numerical prediction in Fig. 1 of Ref. [5j for a 
sphere strictly in the H/2 mid-plane. 



decrease as the confinement becomes stronger. The 
larger positive slopes exhibited by the more needle-like 
spheroids are indicative of motions more sensitive to con- 
finement. Furthermore, the slopes of the same ellipsoids 
in Fig. [6] are larger than those in Fig. [5l indicating that 
diffusion along the ellipsoid short axis is more strongly 
affected by the confinement than diffusion along the long 
axis. Limited comparisons with all available analyti- 
cal and numerical predictions (i.e. the solid curves in 
Figs. [3 [6]), suggest that our data exhibit the gener- 
ally expected trends with increasing confinement. Note 
that solid curves of analytical and numerical predictions 
in Figs, m [6] are for particles forced in the z — H/2 
mid-plane. In real experiment, the measured drag is an 
average at different z. In our experiments, there are 
no detectable interference color changes at the centers 
of ellipsoids. Consequently z-fluctuations are less than 
^ 50 nm ^ H/20. In contrast, numerical results in 
Ref. 5] show that the drag of a sphere at i//3 is very close 
(< 10%) to the drag at H/2. Thus the z-fluctuations of 
our ellipsoids should have negligible effects on particle 
drags. 

Another question that our data holds potential to ex- 
plore concerns the effect of electric double layers on ellip- 
soid diffusion. The electric double layers around charged 
particles in suspension increase their hydrodynamic di- 
ameter and slow down diffusion, especially rotational dif- 
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FIG. 6: (color online) Ratio of theoretical 3D diffusion co- 
efficient [j] along the ellipsoid short axis, D^^ , to the mea- 
sured diffusion coefficient, Dt, for ellipsoids confined at 2b/ H. 
Symbols are the same as those in Fig. [5] Solid curves from 
left to right: theoretical weak confinement predictions for 
aspect ratios a/h — 3, 2, and numerical result Q for aspect ra- 
tion a/b = 1 ranging over both weak and strong confinement 
regimes. 
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FIG. 7: (color online) Ratio of theoretical rotational diffusion 
coefficient 0, -Dg^ to the measured diffusion coefficient, 
De, for ellipsoids confined at 2b/H. Symbols are the same as 
those in Fig. [S] 



fusion because rotational drag is proportional to the vol- 
ume rather than the length of the ellipsoid. This effect 
in rotational diffusion has been observed with depolar- 
ized dynamic light scattering in the regime where ionic 
concentration was low and spheroids small [l9j . In our 
systems such effects are expected to be small due to the 
high ionic strength of the suspension. As can be seen in 
Figs. [U [H and [3 diffusion coefficients are indistinguish- 
able for the samples with 2 mM added salt and no added 
salt. The 6.9 nm screening layer of 2 mM solution lowers 
3D diffusion coefficients by less than 2%. 

Finally, Figure [5] shows the impact of aspect ratio on 
the ratio Da/Dt- Here it is evident that diffusion in 



FIG. 8: (A) Ratio of diffusion coefficients Da/Dfj versus as- 
pect ratio 4> ~ Seven of eight data points were taken 
for 2b/H ~ 0.8; the data point at a/b — 8 was taken with 
2b/ H ~ 0.4. Solid curve: Theoretical result for Da/ Di, versus 
aspect ratio phi =a/b in 3D. (B) Microscope images of the 
corresponding eight particles as a function of aspect ratio (in- 
creasing from left to right). The dimensions of each particle 
from left to right: [2a, 26] = [2.48, 1.2], [2.96, 1.04], [2.8, 0.96], 
[2.88,0.9], [2.64,0.8], [3.12,0.88], [3.3,0.635], [2.4,0.3] fim. 



quasi-2D is quite different from diffusion in 3D. For 3D, 
Da/Db asympotes to 2 at large aspect ratio, as shown by 
the solid theoretical curve. For quasi-2D, Da/D;,, on the 
other hand, grows very rapidly with increasing aspect ra- 
tio Since we expect the stick boundary condition to hold 
in this system, the observation that Da/Dh = 7h/7a > 2 
should be purely due to confinement. A schematic to 
qualitatively capture this basic effect is given in Fig. O 
Imagine the fluid flowing past the ellipsoid. In 3D, the 
fluid flow pathways will be displaced by distances of order 
2b in order to 'go around' the ellipsoid. This fluid flow 
displacement is the approximately the same, whether the 
spheroid is orientated either parallel or perpendicular to 
the flow, and therefore ja and 75 are comparable. In 2D, 
however, the fluid flow pathway displacement is approx- 
imately 26 (or 2a) when the spheroid is oriented parallel 
(perpendicular) to the flow, so that jt/ja diverges with 
a/b. This qualitative picture also explains our observa- 
tion that diffusion along the ellipsoid short axes is more 
strongly affected by the conflnement than diffusion along 
the long axis. Finally, we note that in quasi-2D conflne- 
ment. Da/ Db increases with increasing aspect ratio and 
should eventually saturate Q at a value much larger 
than two, because some fluid will 'leak' between the par- 
ticle and walls. 

In summary, we have found that the anisotropic drag 
coefHcients for ellipsoid diffusion substantially increase 
when the ellipsoids are strongly conflned, especially along 
the short axes. In the future many questions remain 
about these systems will be exciting to explore, includ- 
ing the effects of neighboring ellipsoids and the effects 
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FIG. 9: Schematic of fluid flow around an ellipsoid perpen- 
dicular to the flow. For the drag along short axes, fluid can 
flow around 2b in 3D, but has to flow around long axis 2a in 
2D. 

of other confinement geometries. For example, quasi-lD 
confinement of an ellipsoid will align the ellipsoid along 
the diffusion direction. This effect may compensate the 
drag from boundaries and lead to an optimal diameter 
for ellipsoid diffusion speed in a quasi-lD cylinder. 
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